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Abstract 

If the Higgs boson were the only particle within the LHC accessible range, precision measurement 
of the Higgs's properties would play a unique role in studying electroweak symmetry breaking as 
well as possible new physics. We try to use low energy experiments such as rare B decay to 
constrain a challenging decay mode of Higgs, in which a Higgs decays to a pair of light (~ 1 ~ 2 
GeV) SM singlet S and becomes invisible. By using the current experimental bound of rare decay 
B — )■ Kvv and computing the contribution of i? — t- KSS to (the) B ^ K +0, we obtain an upper 
bound on the Higgs coupling to such light singlet. It is interesting that the partial width of the 
invisible decay mode /i — )• SS hy taking the upper bound value of coupling is at a comparable 
level with h — )■ WW/ZZ or l^/'M^^*) decay modes, making the Higgs identifiable but with a different 
predicted decay BR from [the] standard model Higgs decay. It will then have an impact on precision 
measurement of the Higgs's properties. We also study the implication for cosmology from such a 
light singlet and propose a solution to the potential problem. 
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I. INTRODUCTION 



Searching the Higgs boson, the last missing piece in the Standard Model (SM) of parti- 
cle physics, is one of the essential goals of the CERN Large Hadron Collider (LHC). The 
Minimal Higgs boson model is the simplest solution to electroweak symmetry breaking and 
also the most economic one to be consistent with existing precision measurements. However, 
theoretical considerations suggest that the Minimal Higgs boson model may not be com- 
plete. Being a fundamental scalar, Higgs boson receives quantum corrections of quadratic 
divergence. To solve this, there have been many theoretical proposals which predict various 
new physics at 0{TeV). Direct evidences of new resonances at the LHC can determine what 
the new physics model is. However, if the Higgs boson were the only particle at the LHC 
accessible range, we will have to rely on precision measurements. The precision measure- 
ment of Higgs boson properties can play an important role to confirm electroweak symmetry 
breaking mechanism [1] and test new physics For instance, measurement on top quark 
Yukawa coupling yt is crucial to probe the origin of fermion mass generation while gg ^ h 
production due to the top quark loop directly depends on the coupling yt. On the other 
hand, the contribution from new physics may also change the gg ^ h production rate signif- 
icantly. One interesting scenario will be that at the LHC one does discover the conventional 
Higgs search channels, confirm it is the Higgs and measure its mass but the observed event 
number is much smaller than what we expect for the SM Higgs of that mass. 

However, a new decay mode of Higgs boson that cannot be easily identified will lead to 
the same consequence when the new decay width is comparable with the conventional SM 
Higgs width at the same mass [3]. For instance, if Higgs decay has an invisible mode, it 
is impossible to fu.lly reconstruct such resonance and is very challenging to identify at the 
hadron colliders j^J. 

In this paper, we want to consider the invisible decay of Higgs to a pair of hidden sector 
scalar (S) particles in the minimal extension of the SM 5-^. As the scalar particle is a 
singlet of the SM interactions it can only directly couple to the Higgs by the interaction 
Lagrangian 

—H^HS^ = -H^HS\ (1) 
Ivq 2 

where A is a dimension one coupling constant and vq the vacuum expectation value (vev) 
of the Higgs boson. It is a challenge to identify such invisible Higgs at collider experiments 



2 



TABLE I: Expected BRs in the SM and experimental bounds (90% C.L.) in units of 10^^. The 
SM values for K, vr, K* include the long distance contributions through intermediate on-shell r, 
which can be dominant for it case 11 1. 



mode 


BRs in the SM [9-12] 


Experimental bounds 


B Kvv 


5.1 ±0.8 


< 14 [13j 


B — )• TTUV 


9.7 ±2.1 


< 100 [14] 


B K*uv 


8.4 ± 1.4 


< 80 [15] 


B — )• pvu 


49+0-61 


< 150 [13] 



and obtain any bound on invisible Higgs. The only controlled experiments at this moment 
that can put constraints on such decay mode are through low energy processes such as rare 
B OT K decays. In these processes the Higgs is virtual, not interacting directly to B or K, 
but to top quark and S. Therefore, the only difference is CKM factor, for K it is about 10~^ 
smaller than B, so we would need more than 10^° K^s. Therefore, we just focus on rare B 
decays in this work. 

In Table [U we show the theoretical estimates of branching ratios (BRs) within the SM 
{912] and their current experimental bounds at B factories 13 15] for the decays B 
Mvv. The errors of the SM estimates in Table |T] are mainly due to the hadronic transition 
form factors and the CKM matrix elements, since those decay channels are among the 
cleanest S M p rocesses due to only involving electroweak penguin diagrams [igI, except for 
B — )■ vrz/z/ 111]. Please note that by taking the ratios such as Br(i? — )■ 7rz/z/)/Br(i? — t- vr/z/). 



Br(i? — )■ K*^^) /'Bi[B — >■ plv), we can reduce considerably the uncertainties related to the 
hadronic form factors jl3|- For B — ^ Kvv, similarly one may consider the ratio Br(i? — )■ 
Kvv) I'BiiB — 7- K£^£~) where the uncertainties from the hadronic form factors are canceled 

n 

to a large extent |12| . 

Here we will focus on B^ K^vv decay as its experimental upper bound is closest to 
the SM prediction as shown in Table HI Using the SM expectation value 



BrsM(5+ ^ K^uv) = 5.1 ± 0.8 x 10" 



(2) 



and the current upper bound from BELLE 13| on this final state as 



Br(S ^ K+X) <14:X IQ-^ , 
3 



(3) 



we can derive the corresponding constraint on Higgs invisible decay width. 

To be kinematically allowed in 5 — )■ KSS, the singlet scalar cannot be heavier than 1 — 2 
GeV. Therefore, the scalar can be easily thermalized through the Higgs interactions in the 
early universe. We first discuss its cosmological bound in next section. The third section is 
the discussion on B decay. After taking into all the constraints, we discuss its implication 
in Higgs in the Section IV and finally the conclusion in Section V. 



II. COSMOLOGICAL BOUND AND DECAY OF A HIDDEN SECTOR SCALAR 



If we assume the renormalizability of the theory and allow the mass term quartic self- 
interaction term and the quartic interaction term with the Higgs, the Lagrangian of the 
scalar sector is written as 



^.caiar = ^{dSy - -m'sS' - - ^S'H^H. (4) 

The Lagrangian respects the Z2 symmetry {S — t- —S) thus 5* is a stable particle. Indeed 
this scalar particle can be a good candidate of dark matter. The scalar particle could be 
in thermal equilibrium with the SM sector through interaction with Higgs boson in early 
universe and finally its relic still may survive in the current universe in the form of dark 
matter js], [3]. The relic density is determined by annihilation cross section of the scalar 



particle to the SM particles as 



19| 



^sh ~ 7 r, (5) 



where as is the annihilation cross section of S to the standard model particles through s- 
channel Higgs exchange diagrams and v is relative velocity between annihilating Ss. Since 
we are mainly interested in GeV scale particle, available channels are mainly to light leptons 
(e, /i, (r)) and quarks {u, d, s, c(, b)) and the cross section is obtained as 

{(TsVrel) = r^{ms). (6) 

The precise value of (^(ms) — ~ x'jy^'^ where xj = rrif/ms depends on the actual 

mass of scalar particle and the kinematically allowed channels. We found a stringent con- 
straints on the annihilation cross section considering the WMAP data Vt^h? = 0.1131±0.0034 
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FIG. 1: Cosmological constraints for a stable S from the relic abundance. Allowed parameter space 
in (A = X/vo,ms) plane with rrih = 115, 150 and 185 GeV, respectively. 
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(7) 



If rrih — 150(115) GeV and ms — 1 GeV we get A > 3.5(1.2), respectively, which is within 
the strong coupling regime where the perturbative description of the model is not available. 

In Fig. [H we presented the allowed parameter space in (A = X/vQ,ms) plane by the 
5 year WMAP data on CDM component with various values of Higgs mass (115, 150, 185) 
GeV taking threshold effects into account. Basically a GeV scale mass range, only in which 
range B — )■ KSS is allowed, is not compatible with the cosmological observations ^. On the 
other hand, if the scalar is heavier [ms > 2 GeV) even though the scalar cannot contribute 
to the 5-decays but can be a successful dark matter candidate, if the A coupling is properly 



In Ref. 



2l|, a scalar field in the mass range of 1 GeV has been considered and the authors reached 



the same conclusion with ours: a large coupling constant is required in order to avoid overabundance. 
However, this large coupling constant is ruled out by the B — Ki^P data, as we consider in Sec. III. 
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chosen. 

However, we can easily avoid this cosmological constraint provided that the singlet ac- 
tually decays into light particles since only (absolutely) stable particles can significantly 
contribute to the dark matter density of the current universe. As the longevity of the scalar 
particle is inherited by the Z2 symmetry, a mechanism of breaking symmetry leads to a 
natural way out. Indeed there is a very promising source of the symmetry breaking. Quan- 
tum gravity effect actually allows higher order operators and some of them might break 
global symmetries such as Z2. For instance, the scalar particle may decay to a pair of 
photons or gluons through dimension five operators: 

+ (8) 

where C\ ~ ~ 0(1) are (unknown) parameters. One should notice that both operators 
respect gauge symmetry but break symmetry. The life time of the scalar is suppressed 
by a large cutoff scale (A ~ Mpianck) but certainly much shorter than the age of universe so 
that we can avoid the strong constraint from the relic density measurements. 

III. B KSS AND INVISIBLE HIGGS 

In this section We study the constraint on the interaction term between the Higgs boson 
and the SM singlet from B decays. Specifically we will look at B ^ KSS decay which 
currently has the most stringent experimental upper bound 14 x 10~^ [l^. 

The effective Hamiltonian for this decay can be expressed as 



Hes = —^CJ{1 + l5)bSS . (9) 



Intuitively, b — )■ sSS decay can be divided into two processes: first b quark decays to s quark 
plus a off-shell Higgs boson h, and subsequently h decays to two light singlets. From the 
interaction Lagrangian term \H+HSy2vo, with H+ = (0-, (vq + h - i(tP)/V2), it is easy 
to show that the Higgs boson decay h — t- SS can proceed through a trilinear term XhSS/2. 
But as we will see later, another term Xcj)^ (f)~ S^ / 2vq is also crucial to guarantee the gauge 
independence of the decay amplitude. 

To evaluate the decay amplitude, the Wilson coefficient Cg at scale /i^ = C(mfe) should 
be known, which can be obtained by matching the full theory to the effective theory at 
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scale around mw to obtain Cs{mw) and then evolving down to fib- As the above opera- 
tor does not mix with other effective operators, the QCD running effects can be obtained 



straightforwardly with the calculation of the anomalous dimension of s(l + 75)6 22| : 



= C,{mw) ■ (10) 

\as[mw) J 

Cs{fnw) can be obtained by calculating the diagrams in Fig. 2. Notice that the Higgs boson 
does not couple to s-quark by taking = 0. 

In Fig. 2, the first eight diagrams represent exactly the intuitive picture that first h — )■ sh, 
and then h — )■ SS. Since the later one is a trivial tree level process, one may first focus on 
the construction of an one-loop effective bsh vertex 

= Cbshv:,Vts-s{'^ + i^)hh (11) 



with the coefficient in t'Hooft-Feynman gauge as 23|, |24| 



Cbsh{mw) = TTT^^ — 3 + Xh tz , 12 

where Xt = mf/mly, Xh = m\/m'^ with the approximation m^/ {m^r.m'l^m'f^ ~ 0. Notice 
that this expression is gauge-dependent as the Higgs boson is off-shell. Although the cal- 
culation itself is straightforward, the issues about gauge dependence and renormalization 
scheme ambiguities are a bit subtle which were finally settle down by several groups a few 



years later 251 ]. 



But for the decay amplitude h — > sSS to be gauge invariant, the last diagram in Fig. 
2, i.e. Fig. 2(i), has to be included which (surprisingly at first look) does not contain 
virtual Higgs boson exchange at all. Actually Fig. 2(i) arises from the interaction term 
A0+0~S'^/2t>o. Therefore strictly speaking, h — )■ sSS can not be factorized into b sh and 
h SS. 

Finally, summing all the diagrams, we obtain ^ 



^ This expression has been obtained in 7]- However, in the derivation, they divided the process b — >■ sSS 
into b — > sh and h SS. They then evaluate the bsh vertex with the approximation of vanishing 
Higgs boson mass. But even with these unrigorousness or approximations, they do obtain finally the 
correct expression which is due to the almost completely cancelation between Fig. 2(h) and (i) up to 
0{m1/m\). However, generally this kind of cancelation does not happen and the summation should be 
of order m^/m^j/, as pointed out by Botella and Lim in [25,] . 
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(a) 



(d) 



(g) 




(b) 



'.0 



(e) 



(c) 



(f) 



(i) 





FIG. 2: 6-quark decays to s-quark plus two light singlets. The internal quark lines represent up, 
charm or top quarks, while the internal dashed lines denote Higgs boson (h) or unphysical charged 
goldstone bosons (0). 



The calculation details can be found in the appendix. Here mf, should be evaluated at the 
scale m\Yi but interestingly when combined with the QCD evolution effect of Eq. ffTOl) . one 
has^ 



rrihimw) 



a 



12/23 



mb{mb) . 



(14) 



Please also note that in 10|] authors considered b sSS in an effective theory approach, 
however, with Cg as model independent free parameters. 

To get the decay amplitude, the hadronic matrix element {K~\s{l + 75)6|_B~) is needed 
as input, which can be related to the known form factors through equation of motion. 



{K-\sil + ^,)b\B-) = ^{K-\s^^b\B-) 

nib 

nib V T 



(15) 



We thank the referee for pointing this out to us 
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with the hght-cone sum rules (LCSR) estimation [26| 



2, 0.162 0.173 



^ 1 - g75.4l2 + (1 - §2/5.412)2 
, 0.33 ^^^^ 

^ l-g2/37.46 

n 

As discussed in [26'] , the uncertainty of the dependence of the form factors have not been 
fully analyzed in LCSR but likely to be smaller than that at = which is about 12%. 
Thus as an rough error estimation we assign a universal 12% uncertainty to the above form 
factors. 

Then, the branching ratio can then be obtained 

(niB-niK)'^ 



4m| 



{K-\-s{l + ^,)h\B-)\lq- - 4m| Jl^^^^^^^ - 4m|,. 



and with the values 



27| 



(17) 



Taking as illustration 

mh = 130 GeV, mg = 1 GeV, (18) 



mhirrib) = 4.2 GeV, = 171.3 GeV, A = 0.814 , Ackm = 0.2257 (19) 

and Vts = — AA^KM' "^^^ obtain the branching ratio 

Br{B ^ KSS) = (0.82 ± 0.20) x ( ' ( ' x 10"^° , (20) 

VI GeV/ \ nih J 

where only the form factor uncertainty has been included in the error estimation. 



IV. INVISIBLE HIGGS 

If there exists such light SM singlet scalar, the Higgs decay can be significantly modified. 
For nis = 1 GeV, we take the upper bound on A derived from the B — )■ as 



Br(5 ^ KSS) = 0.82 x 10 



-10 



A \Vl30GeVy 



1 GeV / \ nih J 
< Bre^p(5 ^ K^) - BisuiB Kuu) ~ 8 x lO"'' (21) 
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FIG. 3: Higgs boson decay BR with Invisible decay mode predicted from current upper bound of 
B —7- Ki'u in solid lines. (For comparison, dashed lines are for SM Higgs decay BR.) 

and compute the upper bound of partial width for h — )■ SS. The partial width of Higgs 
decaying into two scalar is 

where A is the dimension one coupling and rrih, ms are the Higgs boson mass and hidden 
sector scalar mass respectively. To illustrate the feature, we scan rrih and plot in Fig. |3] how 
the Higgs decay BR will be changed due to the h — )■ SS decay. The partial width of /i — )■ SS 
is obtained by taking ms = 1 GeV and the A upper bound value computed for that rrih 
point. 

If rrih < 150 GeV, h — )• SS completely dominates the Higgs decay and Higgs is only 
invisible. Even though the traditional invisible Higgs search can be applied to search for 
such modes, it is impossible to identify the resonance through invisible modes at the LHC 

When ruh > 150 GeV, the partial width oi h SS is comparable to the partial widths 
of conventional channels, such as /i — )■ W^W~ oi h ZZ. The multi-lepton searches for 
Higgs resonance are still valid but the decay BRs significantly decrease. If the measured 
event numbers of /i — )■ W^W~ or /i — > ZZ are below the expected numbers. There are 
several possibilities: 
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FIG. 4: Solid lines correspond to the Higgs BR with invisible decay mode predicted from the upper 
bound value for ms = 1 GeV and Br(i? — )■ KSS) = 1 x 10~^ or 1 x 10~^ respectively. Dashed 
lines are the standard SM Higgs decay BR. 

• There are more than one Higgs boson responsible for the W gauge boson mass Mw- 
The vacuum expectation value for the lightest Higgs boson is much smaller than vo so 
that the coupling W^W~H is gv' . 

• The production of Higgs boson is suppressed due to new physics. For instance, gg ^ H 
production is less due to the top quark partner in the triangle loop and significantly 
cancel the top quark loop. 

• There exists unknown Higgs decay mode which cannot be easily identified. Invisible 
Higgs mode that we discuss here falls into this category. Another example is the 



h — )■ uN decay in some TeV neutrino models p]. 

We expect the SuperB or SuperBelle will improve the measurement significantly and 
reduce the allowed region of Br (5 KSS) = Brcxp(-B -> KX) - BrsM(-B -> Kuu). In Fig. 
m we plot how the Higgs decay BR will change accordingly for ms = 1 GeV and improved 
bound on Br (5 — i- KSS). As can be seen, if the value of Br(i? KSS) becomes smaller 
than 2 x 10~^, it will only change the Higgs decay before on-shell WW threshold and won't 
significantly change the heavy Higgs decay. 
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V. CONCLUSION 



We have studied the contribution of virtual Higgs in S — > K^hy assuming Higgs couphng 
to a hght SM singlet scalar S, B ^ KSS. For M5 = 1 GeV, 

BriB KSS) = (0.82 ± 0.20) x (^) ^ (1^) ' x 10- . 

Given the current experimental bound and subtracting the known SM contribution, 

Brexp(5 ^ K^) - BisuiB Kuv) ~ 8 x 10"^ , 

we obtain an upper bound on the coupling between the Higgs and singlet scalar S. We 
take the upper bound value of this coupling and compute the h SS decay partial width. 
It is interesting that the partial width oi h ^ SS decay is at comparable level when the 
Higgs mass is close to the 1^1^ threshold. Consequently, Higgs may still be discovered via 
the conventional Higgs search channels but with a smaller event number. This will have 
some impact on precision measurement of Higgs property. We expect that the SuperB or 
SuperBelle experiments can improve the B — )■ measurement and put a stringent bound 
on possible invisible Higgs decay. 

We have also studied the possible implication in cosmology from this scalar. It turns 
out that for the interesting region of couplings between h and S, such light scalar may not 
have enough annihilation cross section and will then over close the universe if it is a stable 
particle. We propose a scenario where S is not stable in the cosmological scale but only a 
stable particle in B decay or collider environments. 
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Appendix A: Calculation details on 6 ^ sSS 

In the calculations, we use the t'Hooft-Fcynman gauge ^ = 1. p, I and q denote the 
momentum of 6-quark, s-quark and virtual Higgs boson, respectively. We have taken the 
approximation — Tn\ ~ — and dropped a common factor \/m\ in the following expres- 
sions. We get 

Fig.2(a) ^ =^-s(m + (AD 

with xt = m1/m^. 



_2 / 




(A2) 



(A3) 



x+y<l 

with 

Ai(x,y) ~ (1 -x-y)m^ + (x + y)mj . 
The divergence of Fig. 2(b) can be canceled by that of Fig. 2(g): 

7 + ln47r- / dxhi ^ ^— — ^ 

^ Jo fJ' 

The sum of Figs. 2(b) and (g) then gives (taking the scale = mw) 

Fig.2(b + g) = ^^.-(0(1 + 7.)wf^ ''"'"'''"'r''"f'"'''""' . (A4) 

It is clear that for Figs. 2(a),(b),(g), the internal up and charm quarks contributions are 
suppressed at least by m^ J compared to the virtual top quark contribution and can be 
safely neglected. 

For Fig. 2(c), the internal top quark contribution is 

Fig.2(c), ^ =i^ma + T.)6Ml^ ^-'"'-'-f_--/>'--^' . (A6) 
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But here the internal up and charm quarks contributions are not suppressed, which can be 
obtained from the above expression by taking the hmit Xt and changing the correspond- 
ing CKM factors. We then obtain using the CKM unitarity condition, 

Fig.2(c) = + (- 2^?l.^^-.-j3.,-^l)(..-l) _ (^g) 

The virtual top quark contribution to Fig. 2(d) is 
Fig.2(d), = ~'f! + iMph 




dxdv ( 2 In ^^^"^'^^ + i^ + ^ + yX 



(AT) 



with 

A2{x, y) ~ xml + (1 — x)m'^ . (A8) 

The divergence here can be canceled when the contributions from the internal up and charm 
quarks are included, then we get 

For Fig. 2(e), we have 

(AlO) 

Here the internal up and charm quarks contributions are again negligibly small due to the 
O^ml^ Jm^) suppression. It is easy to show that the contribution of Fig. 2(f) vanishes by 
using the equation of motion s{l)l/= 0. The cancelation between Fig. 2(h) and Fig. 2(i) is 
obvious by approximating the Higgs boson propagator i/{q^ — ml) ~ —i/m\. 
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